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Abstract 

We analyze a gradient flow of closed planar curves minimizing the anisoperi- 
metric ratio. For such a flow the normal velocity is a function of the anisotropic 
curvature and it also depends on the total interfacial energy and enclosed area of 
the curve. In contrast to the gradient flow for the isoperimetric ratio, we show there 
exist initial curves for which the enclosed area is decreasing with respect to time. 
J^ '. We also derive a mixed anisoperimetric inequality for the product of total interfacial 

energies corresponding to different anisotropy functions. Finally, we present several 
computational examples illustrating theoretical results. 
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1 Introduction 



k> ' The goal of this paper is to investigate a geometric flow of closed plane curves T ,t > 0, 

H ■ minimizing the anisoperimetric ratio. We will show that the normal velocity /3 for such a 

- - - geometric flow is a function of the anisotropic curvature, the total interfacial energy and 

enclosed area of an evolved curve, 

P = 6{u)k + J^r, (1) 

where k is the curvature and 6{u) > is a strictly positive coefficient depending on the 
tangent angle z/ at a point x G F*. Here J-p is a nonlocal part of the normal velocity 
depending on the entire shape of the curve F = F* and the term 5{v)k represents the 
anisotropic curvature. In typical situations, the nonlocal part is a function of the enclosed 
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area A and the interfacial energy L^ = j^ads, i.e. J-p = J^{A,L^). As an example one 
can consider 

where L = Li is the length of an evolved closed curve F. It is well known that such a 
flow represents the area preserving geometric evolution of closed embedded plane curves 
investigated by Gage [5]. Among other geometric flows with nonlocal normal velocity 
we mention the curvature driven length preserving flow in which /3 = k — ^ J^ k'^ds 
studied by Ma and Zhu [10] and the inverse curvature driven flow preserving the length 
/3 = —k~^ + ^ studied by Pan and Yang [13]. The isoperimetric ratio gradient flow with 
(3 = k — L/{2A) has been proposed and investigated by Jiang and Zhu [9] for convex 
curves and by the authors in [16] for general Jordan curves evolving in the plane. 

The main goal of this paper is twofold. First we derive the normal velocity /3 corre- 
sponding to the anisoperimetric ratio gradient flow. It turns out that P = kcr — Lcr/{2A) 
where k^ is the anisotropic curvature, i.e. /3 has the form of ([I]). We derive and analyze 
several important properties of such a geometric flow. In contrast to the isoperimetric 
ratio gradient flow (c.f. Jiang and Zhu [9], [16] ), we show that the anisoperimetric ratio 
gradient flow may initially increase the total length and, conversely, decrease the enclosed 
area of evolved curves. In order to verify such striking phenomena, an accurate numerical 
discretization scheme for fine approximation of the geometric flow has to be proposed. 
This is the second principal goal of the paper. We derive a numerical scheme based on 
the method of flowing finite volumes with combination of asymptotically uniform tan- 
gential redistribution of grid points. The idea of a uniform tangential redistribution has 
been proposed by How et al in [8] and further analyzed by Mikula and Sevcovic in [TT] . 
The asymptotically uniform tangential redistribution has been analyzed in [131 [12] • The 
scheme is tested on the area-decrease and length-increase phenomena as well as on various 
other examples of evolution of initial curves having large variations in the curvature 

The paper is organized as follows. In the next section we recall the system of governing 
PDEs describing the evolution of all relevant geometric quantities. In section 3 we recall 
basic properties the anisotropic curvature and Wulff shape. We prove an important dual- 
ity identity between total interfacial energies corresponding to different anisotropies. In 
section 4 we investigate a gradient flow for the anisoperimetric ratio. It turns out that the 
flow of plane minimizing the anisoperimetric ratio has the normal velocity locally depend- 
ing on the anisotropic curvature and nonlocally depending on the total interfacial energy 
and the enclosed area of the evolved curve. Section 5 is devoted to the proof of a mixed 
anisoperimetric inequality for the product of two total interfacial energies corresponding 
to two anisotropy functions. In section 6 we investigate properties of the enclosed area 
for the anisoperimetric gradient flow. In contrast to a gradient flow for the isoperimetric 
ratio, we will show that there are initial convex curves for which the enclosed area is 
strictly decreasing. Finally, in section 7 we construct a counterexample to a comparison 
principle showing that there initial, one of them contained in the closure of other one, 
and such that they intersect immediately when evolved in the normal direction by the 
anisoperimetric ratio gradient flow. Section 8 we derive a numerical scheme for solving 
curvature driven flows with normal velocity depending on no-local terms. The scheme is 
based on a flowing finite volume method combined with a precise scheme for approxima- 
tion of non-local terms. We present several numerical examples illustrating theoretical 
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results and interesting phenomena that are present in the gradient flow for anisoperimetric 
ratio. 

2 System of governing equations and curvature ad- 
justed tangential redistribution 

In this section we recall description and basic properties of geometric evolution of a closed 
plane Jordan curve F which can be parameterized by a smooth function x : [0, 1] — t- M^ 
such that r = Image(a;) = {x{u); u G [0,1]} and \dux\ > 0. We identify the interval 
[0, 1] with the quotient space M/Z by imposing periodic boundary conditions for x{u) at 
u = 0,1. We denote d^F = dF/d^, and \a\ = y/a ■ a where a ■ 6 is the Euclidean inner 
product between vectors a and b. The unit tangent vector is given by T = dux/\dux\ = 
dgX, where s is the arc- length parameter ds = |9ua^|dM. The unit inward normal vector is 
deflned in such a way that det(T, N) = 1. Then the signed curvature k in the direction 
N is given hj k = det{dsX,dgX). Let z/ be a tangent angle, i.e., T = (cos i^, sin z/)""" and 
N = {— sin u, cos u)"^. From the Frenet formulae dgT = kN and dgN = —kT we deduce 
that dgV = k. 

Geometric evolution problem can be formulated as follows: for a given initial curve 
F° = Image(cc°) = F, flnd a family of curve {F*}t>o, F* = {x{u,t)] u G [0,1]} starting 
from F° and evolving in the normal direction with the velocity /3. In this paper we follow 
the so-called direct approach in which evolution of the position vector x = x{u,t) is 
governed by the equation: 

dtX = (3N + aT, x{-,0) = x%-). (2) 

Here a is the tangential component of the velocity vector. Note that a has no effect on the 
shape of evolving closed curves, and the shape is determined by the value of the normal 
velocity (3 only. Therefore, one can take take a = when analyzing analytical properties 
of the geometric flow driven by (|2]). On the other hand, the impact of a suitable choice of 
a tangential velocity a on construction of robust and stable numerical schemes has been 
pointed out by many authors (see [13 [IS] and references therein). 

In what follows, we shall assume that /3 = 5{v)k + J-r where 5{v) > is a strictly 
positive 27r-periodic smooth function of the tangent angle v and J-p is a nonlocal part 
of the normal velocity depending on the entire shape of the curve F. According to [12] 
(see also [TT|[T3]) the system of PDEs governing evolution of plane curves evolving in the 
normal and tangential directions with velocities /3 and a reads as follows: 

dtk = dlP + adsk + A;2/3, (3) 

dtu = dsP + ak, (4) 

dtg = {-k/3 + dsa)g, (5) 

dtx = 5{u)d'^^x + adsx + JrAT, (6) 

for u G [0, 1] and t > 0. Here g = \dux\ is the so-called local length (c.f. [H]). A solution 
to ([3])-([6]) is subject to periodic boundary conditions for g,k,xatu = 0, 1, z/(0, t) = z^(l, t) 
mod(27r) and the initial condition k{-,0) = ko{-),h'{-,0) = i^o{-),g{-,0) = go{-),x{-,0) = 
cco(-) corresponding to the initial curve F° = Image(a;°). 
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Local existence and continuation of a classical smooth solution to system ([3])-([6]) has 
been investigated by the authors in [T5l [T6] . In this paper we therefore take for granted 
that classical solutions to ©-(E]) exists on some maximal time interval [0,Tmax) (cf. 
[I61[T3]). 



3 The WulfF shape and interfacial energy functional 

The anisotropic curvature driven flow of embedded closed plane curves is associated with 
the so-called interfacial energy density (anisotropy) function a defined on F. It is assumed 
that a = a{u) is a strictly positive function depending on the tangent angle u only. With 
this notation we can introduce the total interfacial energy 



L,(r)= f a{u)ds 



associated with a given anisotropy density function a. If a = 1 then i^i(r) is just the total 
length L[r) of a curve F. The Wulff shape is defined as an intersection of hyperplanes: 

Wo- = ll {x = {xi,X2)^; —x-N = xismiy — X2Cosiy<a{u)^. 

If the boundary dW^ of the Wulff shape is smooth and it is parameterized by dW„ = 
{x = —a{v)N + a{i')T, v G [0, 27r]}, then, it follows from the relation d^v = k that 

T = dsX = {-a'{u) + a{iy))kN + (a(z/) + a\v))kT. 

Hence a{u) = cr'{iy) and (cr(i^) + a"{u))k = 1 holds and the boundary dWa can be param- 
eterized as follows: 

dW^ = {x; x = -a{u)N + a\u)T, v E [0, 27r]} , 

and its curvature is given by k = (cr(i^) -|- o-"{iy))^^. Let us denote by k^r the anisotropic 
curvature defined by k^ '■= {o'{i') + a'\v))k. It means that the anisotropic curvature k„ 
of the boundary dW^ of the Wulff shape W^ is constant, k^ = I- Moreover, the area 
\Wa\ = A{dWa) of the Wulff shape satisfies: 



2 JdW„ 2 Jqw„ 2 



^ Jdw„ ^ Jaw, 

Clearly, \Wi\ = vr for the case a = 1. If we consider the anisotropy density function 
a{i') = 1 + e cos(mi/) for m = 2, 3, ■ • • , £(m^ — 1) < 1 then the area of W^ can be easily 
calculated: 

\W^ = \j a{u)ds = \r a{a" + a)du=l{2-e\m'-l)). (7) 

In Figure [1] we plot shapes of dW^ for various degrees m. 

Since the global quantities evaluated over the closed curve F do not depend on the 
tangential velocity a we may take a = 0. Hence dtg = —kf3g and dti' = dgfS. These 
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m = 2 



m, = 3 



m = A 



m = 5 



7n = 6 



Figure 1: The Wulff shapes W^ for m = 2, 3, 4, 5, 6 and e = 0.99/ {m"^ - 1). 



identities follow from (j4j) and ([5]) with a = 0. Recall that dsiy = k. Therefore ds<j'{i') = 
a"{h')dsiy = a"{u)k and so j-^a"{u)kds = 0. Hence 

/Co-ds = / akds 
r ^r 

holds. For the time derivative of Jp k^ds we obtain 

— / k^ds = -77 / (^kgdu = / [dt{(Tk)g + akdtg\du = / [dt{(Tk) - (Tk'^f3]ds 

[kdtaiu) + a{iy)dtk - a{u)k^/3]ds = f[ka'iu)dtu + (j{v){dtk - P/3)]ds 
[ka\v)dsP + a{v)dlP]ds = 0, 

because dtk — k'^f3 = dg/3 and ka\u) = a\u)dsi' = dga^u). From the previous equality we 
can deduce the following identity: 



krrds 



r* 



kads, for any < t < T„ 



ro 



where the family of planar embedded closed curves F*, t G [0, T^ax), evolves in the normal 
direction with the velocity /3. 

Now, let us consider an evolving family of plane embedded closed curves F*, t G [0, T], 
homotopicaly connecting a given curve F = F° and the boundary F-^ = dWa of the 
Wulff shape W^- The homotopy can be realized by taking a suitable normal velocity /3 
(eventually depending on the position vector cc). Using such a normal velocity we deduce 
the identity: 



fk„ds= I kAs = L{dW,). 

Jv JdWa 



(9) 



It means that Jp ka-ds is equal to the length of the boundary dW„ of the Wulff shape. 
The same result has been recently obtained by Barrett et al. in [21 Lemma 2.1]. We 
can say that identity ([9]) is a generalization of the rotation number: -^ Jp kds = 1, since 
2n = L{dWi). 
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Remark 1 Identity ^ can be easily shown for convex curves. Indeed, ifV is convex then 
its arc-length parameterization s can he reparameterized by the tangent angle v G [0,27r]. 
We have dsiy = k > and therefore ds = k^^du. Hence 

kads = / akds = / cr(i^)dz/. 
'r ir Jo 

For the length L{dW^) of the boundary of a convex Wulff shape we obtain 
L{dW^) = ds= T^^= I [<y{^) + (r"{i^)]dv = I a{v)dv. 

JdWa Jo ^ Jo Jo 

Therefore f-pk^ds = L{dWa) because /Q'^cr"(i/)dz/ = and k = [cr(z/) + a'^u)]"^ on dW^- 
If r is not convex we can apply the famous Grayson's theorem ^. We let it evolve 
according to the normal velocity P = k until a time t = T when T^ becomes convex. Using 
^ and previous argument we again obtain identity ^. 

Let us denote by Li the total interfacial energy corresponding to a = 1, i.e. Li = L. 
Let r = dWi be the unit circle. Then, by applying identity ([9]), we deduce 

Li{dW^) = L^{dWi). (10) 

Latter identity can be rephrased as follows: the length of the boundary dW^ of the Wulff 
shape equals to the total interfacial energy of the unit circle. It can be easily generalized 
to the case of arbitrary two anisotropies (j(z/) and /i(z/). We have the following proposition: 

Theorem 1 Let a and fi be two smooth anisotropy functions satisfying cr(z/) + ^"{u) > 
0,/i(z/) + fi"{u) > 0. Then the duality 

L^idW^) = L^idW^) (11) 

between total interfacial energies of boundaries dW^ and dW^ of Wulff shapes holds. 

Proof. Notice that the Wulff shapes W„ and W^ are convex sets because a{u)+a"{u) > 
and ^{v) + ^"{v) > hold. For the curvature k at the boundary dW^j we have k = 
[a{p) + a"{u)]-^ and so 

p nliv 1 /'27r 

L^{dW,)= fiiiy)ds= /i(z/)-dz/= / fi{u){aiiy)+a"{iy))du 

JdWa Jo ^ Jo 

f.27r 

[/i(z/)(T(i.) - a'iiy)fi'{iy)]du = L^idW^), (12) 

'0 

arguing vice versa. <) 



On a gradient Bow of plane curves minimizing the anisoperimetric ratio 1 

4 Gradient flow for the anisoperimetric ratio. 

Recall that for the enclosed area A = A(r^) and the total length L = L(r*) for a flow of 
embedded closed plane curves driven in normal direction by the velocity /3 we have 

— A+ / f3ds = 0, :r^+ / k/3ds = 0, (13) 

dt Jpt at J-pt 

(c.f. [H]). Using governing equations (I3l)-(l6]), for the total interfacial energy L^- = L„{T^) 
of a curve F*, we obtain 

—L^ = — / a{u)ds = -Tj. cr{u)gdu = / [a'{u)dtug + a{u)dtg]du 

[cT'{u)d,f3 - a{iy)kf3]ds = - f[a"{u)dsiyp + a{u)kf3]ds (14) 

= - /[^"(^) + cr(z/)]A;/3ds = - f KPds. 

Here we have used the governing equations ([5]) and (j4j) (with a = 0) and the identity 
dsV = k. For the anisoperimetric ratio 

n,(r)- ^"W 



wMi^y 



we have no-(F) > 1 and, in particular, Ilu{dWa) = 1 (see Remark[2]). Taking into account 
identities (ITU) and flT3l) we obtain 



d_^ L^dtL, LldtA L 



TaIX''-T^'''- 



dt " 2\WM ^\W,\A^ 2\W, 

Hence, the flow driven in the normal direction by the non-locally dependent velocity 

/3 = ^. - ^ (15) 

represents a gradient flow for the anisoperimetric ratio Hg- with the property djla < for 
/3 ^ 0. Notice that /3 = on F if and only if F oc dWc,-, i.e. F is homotheticaly similar to 

In the case cr = 1 the isoperimetric ratio gradient flow has been analyzed by Jiang 
and Zhu in ^ and by the authors in [T3]. In this case the normal velocity has the form: 

I3 = k-L/{2A). 

5 A mixed anisoperimetric inequality 

The aim of this section is to prove a mixed anisoperimetric inequality of the form 



A{V) 



> i^.,M, (16) 



which holds for any C^ smooth Jordan curve F in the plane. Here K„^^ > is a constant 
depending only on the anisotropy functions a and fj, such that cr{h') + a"{iy) > and 
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/i(i^) + /^"(i^) > hold for any u. The idea of the proof is rather simple and consists in 
solving the constrained minimization problem: 

mmL,(r), s.t. L^{T) = cA{T), (17) 

where c > is a given constant. To this end, let us assume that a curve F = r(a;) is 
parameterized by a C^ smooth function x : S*^ — t- M^. If we denote g = g{x) = \dux\ 
the local length then, for the derivative of g in the direction y : 5*^ — )■ M^, we obtain 
2g{x)g'{x)y = 2{duX ■ duV) and so 

g\x)y = {T-dsy)g. (18) 

Here and here after, for scalar-valued function f{x) and vector-valued function f{x) = 
{fi{x), f2{x))'^ we denote their derivatives in the direction y by 

f {x)y := Vf{x) . y = hm , / {x)y := ^ j^j^ 

respectively. 

As for the tangent vector T = T(x) = (cos z/, sin z/)""" we have T{x) = g~^duX and 
so T'{x)y = g-'duy - g~'duXg'{x)y = d^y - (T ■ dsy)T = {N ■ dsy)N. As N = 
(— sinz/, cosz/)^, for the derivative of the tangent angle u = I'^x), we obtain 

iy'{x)y = N ■ d,y. (19) 

Recall that k^ := (o-(z^) + a'\u))k and ds^ = k. Since Lfj{T) = f-pCrds = f^ a{iy)gdu we 
obtain 

K{nx))y = j [a\u)u\x)yg + a{u)g\x)y] du = j [a'{u){N ■ d,y) + a{u){T ■ d^y)] ds 

= - j [a'\u)dsu{N ■ y) - a\u)k{T ■ y) + a\u)dsu{T ■ y) + a{u)k{N ■ y)] ds 

= - j [a{u) + a"{u)] k{N ■ y)ds = - j k^{N ■ y)ds. 

Hence 

LUr{x))y = - j k^{N ■ y)ds, L'^{V{x))y =- j k,{N ■ y)ds. (20) 

For the area A = A{T) enclosed by a Jordan curve F = T{x) we have A = ^ f^ det(a;, dux)d 
Therefore 

A'{T{x))y = - / det{y,dux) + det{x,duy)du= / det{y,dux)du = det{y,T)ds. 
2 Jo Jo Jv 

Since det(j/, T) = —y ■ N we obtain 

A{T{x))y = - [ N-yds. (21) 



u. 



r 



In order to solve the constrained minimization problem (fT71) we introduce the Lagrange 
function C{x, A) = L^iT{x)) + A(L^(F(a;)) - cA{T{x))) with A > 0. 
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Then the first order condition for T = T{x) to be a minimizer of (IT7j) reads as follows: 
= C'^{x, X)y = L'^{T{x))y + \{L'^{T{x))y — cA'{T{x))y) sX x = x. Latter equality has 
to be satisfied for any smooth function y : S^ ^ M?. Taking into account ( l20l) and ( l2Tll 
we obtain 

kfj + Xk^ = Ac, on f . 

It means that 

k^ = Ac, on r, where a = a + Xfi. (22) 

In other words, f = j;:dWa (up to an afiine translation in the plane M^). The Lagrange 
multiplier A G M can be computed from the constraint L^{r) = cA(r). It follows from 
duality (ITT]) (see Proposition [T]) that 

L,{dW,) = L,{dW,) = L^idW,) + XL^idW,) = L^idW,) + 2XA{dW^). 

To calculate the enclosed area A(T) = j^A{dWa) we make use of the identity A{dWa) = 
^L^{dWa). Clearly, as o" = o- + AyU we obtain 

MdW^) = L,{dW-„) + XL^{dW^) = L^idW^) + XL^idW^) 
= L,{dW„) + XL^{dW,) + XL,{dW^) + X^L^{dW^) 
= 2A{dW„) + 2XL„{dW^) + 2X^A{dW^). 

Since ^L^(aiys) = L^(f ) = cAif) = ^A{dW^) we end up with the identity 
^ (L^idW^) + 2XA{dW^)) = -^ {A{dW^) + XL^{dW^) + X''A{dW^)) . 

AC A C 

Since the Lagrange multiplier A > it is given by A = ^y A{dWcr) / A{dW^) . 
Furthermore, 

L^idW^) = L^idW,) = L^idW^) + XL^idW^) = 2A{dW^) + XL^idW^). 

Now, let r be an arbitrary C^ smooth Jordan curve in the plane. Set c = L^(r)/y4(r). 
Then 

L.(r)L^(r) _ ^^ ,^, ^^^ ,-,_ c 



^.^. cL.(r) > cL^(r) = Y^L^idW,) = 2^A{dW^)A{dW,) + L^{dW^). 

In summary, we have shown the following mixed anisoperimetric inequality: 
Theorem 2 Let T be a C"^ smooth Jordan curve in the plane. Then 



^.^. >i^.,M, where K^,, = 2^\W^\\W,\ + L,{dW^). (23) 

The equality in ( fl^) holds if and only if the curve T is homothetically similar to the 
boundary dW^ of a Wulff shape corresponding to the mixed anisotropy function a = 

v^wg^ + VW/i. 

Remark 2 If a = fi = 1 we obtain the well known isoperimetric inequality L(r)^/74(r) > 
Kii = 2^/7r^+L{Wi) = An. If a = n we obtain the anisoperimetric inequality L^{Ty / A{T) > 
Ka,a = 2^J\Wa\'^ + L^{dWfj) = 4:\Wfj\. Finally, if fi = 1 we obtain the mixed anisoperi- 
metric inequality 
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6 Convexity preservation. Temporal area and length 
behavior 

In this section we analyze behavior of the enclosed area A{r^) of a curve F* evolved in the 
normal direction by the anisoperimetric ratio gradient flow, i.e. 13 = k^ — L„/(2A). 

First we prove the preservation of convexity result stating that the anisoperimetric 
ratio gradient flow preserves convexity of evolved curves. In the case of the isoperimetric 
ratio gradient flow of convex curves with f3 = k — L/{2A), the convexity preservation 
has been shown by Jiang and Pan in [9]. However, similarly as Mu and Zhu in [TO] . 
they utilized the Gauss parameterization of the curvature equation ([3]) by the tangent 
angle u and this is why their results are applicable to evolution of convex curves only. 
In our paper we first prove convexity preservation based on the analysis of the curvature 
equation (^ with arc-length parameterization. Moreover, we show the anisoperimetric 
ratio gradient flow may initially increase the total length and decrease the enclosed area. 
This phenomenon cannot be found in the isoperimetric ratio gradient flow (c.f. [9| [15]). 

Theorem 3 Let F*,t G [0,Tmax), be the anisoperimetric ratio gradient flow of smooth 
Jordan curves in the plane evolving in the normal direction by the velocity /3 = /Cg- — ff ■ 
// the curve T^° is convex at some time to ^ [0,Tmax) then F* remains convex for any 

t £ rO) Tmax)- 

Proof. Since dfk = d'^P + k'^P, dtv = dgP = dsk^ and (3 = k„ + J-p we have 
dtk^ = 6{u)dtk + 6'{u)kdtu = 5{v)dlk, + <5(z/)P/3 + 6'{u)kdsk^ 
= 5{v)dlk„ + -^kl(3 + 5\v)kdsk„, 

where 5{v) := cr(z/) + cr"{iy) > 0. 

Let us denote by K(t) = minpt A;(j(., t) the minimum of the anisotropic curvature 
kcr = 6{v)k over the curve F*. Denote by s*{t) G [0, L*] the argument of the minimum of 
k„, i.e. Kit) = k^{s*{t),t). Then dsk^{s*{t),t) = and d'^k^{s*{t),t) > 0. Hence 

K'{t)>^^K{tr{K{t)+T',), 

where A(t) = 6{u{s*{t),t)), and J"^ = -Ll/{2A^). Notice that A(t) > A™„ := min^ (5(z/) > 
for all t G [0, T^ax)- Suppose that i^ is a solution to this ordinary differential inequality 
existing on some interval [to,Tmax) and such that K{to) > 0. Then, it should be obvious 
that K(t) > for t G [to,Tmax) provided that 

inf J"^ > -oo, (24) 

to<t<t* 

for every < t* < Tmax- In order to prove convexity preservation for the anisoperimetric 
ratio gradient flow it is therefore sufficient to verify that the nonlocal part J-'p = — L^/(2A*) 
remains bounded from below for t > to- To prove boundedness of J-'p from below we utilize 
a property of the anisoperimetric ratio. Indeed, as [3 = k^ — Lfj/{2A) represents gradient 
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flow for the anisoperimetric ratio H^ = (L^)^/(4|lVo-|^*) we have 1 < H^ < n[^ for all 
t e [0,Tmax)- Thus 



Now, since the classical solution exists on the time interval [0, Tmax) then info<j<i* L^ > 
for each < t* < Tmax and the estimate ( 12^ follows. (} 



In what follows, we shall investigate the enclosed area and length behavior of curves 
evolved by the normal velocity [i = k^ — Lcr/(2A) representing thus a gradient flow for 
the anisoperimetric ratio. 

Using the area equation (fT3l) we obtain 

±A^-l p,s ^ - I {k - ^) d. ^ -LiSW,) + ^. (25) 

By applying the isoperimetric inequality (see Remark |2]) for the case a = 1 and any curve 
r = r*, the following inequality: 

dt 2A{r) - 

holds. It means that the gradient flow for the isoperimetric ratio does not decrease the 
enclosed area. On the other hand, if anisotropy density function a ^ const, then for a 
curve r = r* oc dW^ corresponding to the Wulff shape W^ with the anisotropy function 
a = a/tt a + ^llVo-l we obtain 

^A = -LidW^) + ^ = -LidW^) + i/C,i = v/^tlg - hidW^) < (26) 

due to the isoperimetric inequality L{dWa-)'^ > iiiA^dWa-) = inlWa-l- It means that the 
gradient flow for the anisoperimetric ratio may initially decrease the enclosed area for 
special initial curves. 

Next we recall the isoperimetric inequality by Gage. According to jl] the following 
inequality holds: 

k^ds > 71- (27) 

r - A 

for any convex C^ smooth Jordan curve in the plane. The equality in fl27|) holds iff F is 
a circle. Therefore, in the case of isoperimetric gradient flow with a = 1 and the convex 
curve r*, we have 

-^L* = - /" A;/3ds = - /" Pds + TT^ < 0. (28) 

However, if a ^ const is a smooth nonconstant anisotropy density function such that 
a + a" > 0, there exists an initial curve F" such that the length L* may initially increase, 
i.e. ^L* > at t = 0. Indeed, let F'' be an initial curve which homothetically similar 
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to the boundary dW^ of the Wulff shape corresponding to the mixed anisotropy function 
a = aa + b where a,b > are constants. Then 

o a a 

Hence k^ = ^ — ^k on the Wulff shape F = r° = dW^ because k^ = 1 on dW^- Using 
( ITS]) , we have 



dt Jt A{dW^) aJr a A{dW^) 

at t = 0. Since a = {a — b)/a we have 

L^idW^) = [ ^^ds = -L^idW^) - -L{dW^) = -A{dW^) - -L{dW^ 
Jqw^ a a a a a 



Thus 



^L^='-(l fcM. - .^i?5^ 1 > att = 0, 



dt a \J9w, A{dW^] 

due to inequahty (1271) and the fact that dW^ is a convex curve different from a circle for 
a ^ const. 

In summary, we have shown the following result. 

Theorem 4 If a = 1 then the isoperimetric ratio gradient flow with the normal velocity 
P = k — L/(2A) is area nondecreasing and length nonincreasing flow of smooth Jordan 
curves r*,t G [0,Tmax) in the plane provided that r° is a convex curve. 

Assume the anisotropy function a is not constant and such that a + a" > 0. Let r° 
be an initial curve which is homothetically similar to the boundary dW^ of a Wulff shape 
with the modified anisotropy density function a = aa + b where a, b are constants, a,b > 0. 
Then, for the anisoperimetric ratio gradient flow r*,t G [0,Tmax), evolving in the normal 
direction by the velocity (3 = k^ — ^, we have 

1. ftL{T^) >0 att = 0. 



2. If, moreover, a/b = \pKl\/\W„\ then ;^^(r*) < at t = 



7 A counterexample to the comparison principle 

A key role in the proof of the famous Gage-Hamilton-Grayson theorem on shrinking 
curves evolved by the normal velocity /3 = /c is played by the comparison argument. It 
states that two smooth curves, one of them included in the closure of the interior of the 
second one, evolved by the normal velocity P = k never intersects each other O [7]. A 
analogous comparison property holds for the so-called afiine scale geometric flow in which 
/5 = /c3 (c.f. PP). The aim of this section is to demonstrate that the analogous comparison 
property does not hold for the anisoperimetric gradient flow. 
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Clearly, any curve F homotheticaly similar to the Wulff shape dW^j is a stationary 
curve, i.e. /3 = on T. Indeed, k^ = 1 on dW^ and J^ow^ = -L^{dWa)/{2A{dW^)) = -1 
and therefore 13 = k^ + J-'aw^ = on dW^- 

In what follows, we shall construct a smooth initial curve F" containing in its interior 
the Wulff shape dW^ and such that F* intersects the stationary Wulff shape dW^ for all 
sufficiently small times < t <^ 1. The construction is as follows. First, we shall construct 
a nonsmooth curve F as the union F = dW„ U r ■ dWi of the Wulff shape dW^ and the 
circle r ■ dWi of a radius r > touching the Wulff shape from outside at a point y (see 
Figured (left)). For such a curve we have 



^f 



L,(f) L^idW^) + L,{r ■ dW,] 



2A(f) 2{A{dW^) + A{r ■ dWi)) 



L,{dW^)+rL,idW,: 
L^{dW^) + 27rr^ 



Lo-(F)/2. Hence J-f > —1 provided that the radius r is sufficiently 



1 + ecos(m.h'), then r > 1 because 



because A{dWa) 

large, r > L„{dWi)/2n. For instance, if (j(z/ 

Lcr{dWi) = 2tt (see section 3). 

Let x° G F n dW„ be a point different from y and belonging to a part of the curve 
representing the Wulff shape. Now, let us construct an initial smooth curve F° which is 
a continuous perturbation of F, it contains the Wulff shape in the closure of its interior, 
and such that F° = F in some neighborhood 0{x) of x (see Figure |2] (right)). The 
anisoperimetric ratio gradient flow F*,t > 0, starting from F" intersects the stationary 
Wulff shape dW^j in the neighborhood 0{x) for any time < t <^ 1 (F* is plotted by 
a dashed curve in Figure |2] (right)). This is a consequence of the fact that the normal 
velocity /3 at x^ is strictly positive for t = because /3 = k^ + J^ro = 1 
and f n 0{x) = dW^ n 0{x), i.e. k„ = I ai x^. 



fo 



J~YQ ^ U cLX X , 



Y = dW^yjr-dW, 



r-dW, 




Figure 2: An initial nonsmooth curve F (left) and its smooth perturbation F° 
Failure of a comparison principle occurs at the point x^ and < t <^ 1. 



fright). 



8 Numerical experiments 



In this section we present several examples of evolution of plane curves minimizing their 
anisoperimetric ratio. Our scheme belongs to a class of boundary tracking methods taking 
into account tangential redistribution. In construction of the scheme we employed flowing 
finite volume discretization in space with a nontrivial tangential velocity, and semi-implicit 
discretization in time. The advantage of applying a nontrivial tangential velocity consists 
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in its capability to overcome various numerical instabilities of a flow of plane curves like 
swallow tails and/or merging of numerical grid points leading to break-up of the numerical 
scheme known for the case when the numerical scheme is constructed with no tangential 
redistribution. One can find recent progress in [151 [IS]- Such a scheme is simple and 
fast, but even if the original problem has a variational structure, it is unclear that the 
discretized problem has variational structure. On the other hand, in [3], the authors 
proposed a semi-discrete scheme with variational structure in a discrete sense. Their 
scheme has the second order of accuracy in time. However, discretized polygonal curves 
are restricted to a certain class of curves which is analogous to the admissible class in 
crystalline curvature flows or crystalline algorithm. In what follows, we propose a hybrid 
scheme taking into account advantages from both aforementioned schemes. 

Discretization scheme. For a given initial A^-sided polygonal curve V^ = IJj=i'^i^' 
we will find a family of A^-sided polygonal curves {V^}j=i^2,---, V^ = [Ji^iSj, where 
Sj = [xl_i,x-l] is the i-th edge with Xq = x\ for j = 0, 1,2, ■ ■■. The initial polygon 
V^ is an approximation of T^ satisfying {a;°}^^ C V^CiT^, and V^ is an approximation of 
r* at the time t = tj, where tj = jr is the j'-th discrete time (j = 0, 1, 2, ■ ■ ■ ) if we use a 
fixed time increment r > 0, or tj = J^iZo '^i (j = 1; 2, ■ ■ ■ ; to = 0) if we use adaptive time 
increments t; > 0, / = 0, ■ ■ ■ , j — 1. The updated curve P-'+^ is determined from the data 
for V^ at the previous time step by using discretization in space and time. Our two steps 
scheme will be constructed as follows: in the first step we construct moving polygonal 
curves which is continuous in time and discrete in space. In the second step we make use 
of the semi-implicit time discretization scheme for moving polygonal curves. 

Step 1: Moving polygonal curves. Let V{t) = [J^^i Si{t) be an A^-sided polygonal 
curve continuously in time with V{0) = V^, where Si{t) = [xi^i(t),Xi(t)] is the i-th edge 
and Xi(t) is the z-th vertex {i = 1,2, ■ ■ ■ , N; Xo(t) = XM(t)). The length of Si is denoted 
by Tj = \xi — Xi_i\. The i-th unit tangent vector Tj can be defined as Tj = (ajj — Xi_i)/ri, 
and the i-th unit inward normal vector Ni = T^^, where (a, b)-^ = {—b, a). Then the i-th 
unit tangent angle Ui is obtained from Tj = (cosz/j, sini/j)""" in the following way: Firstly, 
from Ti = (Til, ^12)""", we obtain 1^1 = — arccos(Tii) if T12 < 0; z/i = arccos(Tii) if T12 > 0. 
Secondly, for i = 1, 2, ■ ■ ■ , A^ we successively compute z/j+i from z/j: 

z/j — arccos(/), ii D < 0, 

z/j + arccos(/), ii D > 0, where D = det(Ti, Tj+i), / = Tj ■ Tj+i. 

z/j, otherwise. 

Finally, we obtain z/q = z/i — (z/^r+i — z/^v). Then the i-th unit inward normal vector A^j is 
Ni = (— sin z/j, cos z/j)^. 

Let us introduce the "dual" volume S* = [x*, Xi] U [xi, x*_^_-^] of Si, where x* = {xi^i + 
a;j)/2 is the mid point of the i-th edge Si {i = 1,2,- ■■, N; x*J^_^_-^ = a;*). The length of 
S* is r* = (rj + rj+i)/2. Then the total length of P is L = X]i=i '^j — 12i=i ^i^ ^^^ ^^^ 
enclosed area oi V is A = — J2i=ii^i ' -^i)'^j/2 = J2i=i ^i~~i ' ^i/2- 

We define the i-th unit tangent angle oiS* by z/* = (z/i + z/i+i)/2 = z/j+0j/2, where 0j = 
z/j+i — z/j is the angle between the adjacent two edges. Then the i-th tangent vector at the 
vertex Xi is T* = (cos z/*, sin i/*)'^ and the inward normal vector A^j* = (— sin u*, cos i^*)^. 
Hereafter we will use the following abbreviations: 

Cj = cos---, Sj = sm--- (z = 1,2, ■ ■■ ,A'). 
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Then it is easy to check that 

T* = cm + SiNi, N* = CiNi - sm (i = 1, 2, ■ ■ ■ , iV). 
The evolution equations of V{t) read as follows: 

Xi = a,T: + (3iN* (z = l,2,---,iV), (29) 

where ctj and /?, are quantities defined on S*. Here and hereafter, we denote it = du/dt. 
The tangential velocities {aj} are defined below and the i-th normal velocity /3j is defined 
such as 

A = ^^^ (^ = l,2,---,iV), (30) 

where the i-th normal velocity (3* is defined on Si. It is an approximation of ([I]) such as 

P: = 6A + J'r (z = l,2,--- ,iV). 
Here ki is the i-th curvature and the constant value on Si defined as 

^__tan(^./2) + tH^-i/2) (; = i, 2. . . . . ^r), 

ri 

which is the same as the polygonal curvature in [3], and 6i is an approximation of Slui) 
defined later. Then we obtain the time evolution of the total length of V{t): 



N N 

i=l i=l 



L = -2j]As^ = -$^M*r 



and the time evolution of the enclosed area of V{t): 



N N N 

2_^ PiCiT, + 2^ aiSi 

i=l 1=1 

N ^ , N 






.,. ^+1 -2ri + ri_i v^ „ ^ n+i - Ti 



Pi T + 2^ «iS* 



These identities represent a discrete version of equations ( 1T3|) provided that the distribu- 
tion Tj = L/A^, 2 = 1, 2, ■ ■ ■ , A^, is uniform because the error term err^ = is vanishing. 
To realize this uniform distribution asymptotically, we assume that 



n - ^ = Vie ^^*^ 



[YliLi Vi = 0, f{t) ^ oo as t ^ T„,ax < oo 
By using a relaxation term uit) = f'(t) we obtain 



f / r \ pTmax 

r,--=(--rAuit), J uit)dt = oo (z = 1, 2, ■ ■ ■ , iV). (32) 

Taking into account the relations: 

fj = {xi - Xi^i) ■ Ti = -/3iSi - A-iSj-i + Citti - Ci-itti-i = J? ~^ [j? ~'^i] ^(^)' 
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we deduce A^ — 1 equations for tangential velocities Oj (z = 2, 3, ■ ■ ■ , A^): 
tti = — H tti (i = 2,3, ■ ■■ , A^), 

^i = ^2 + ^^3 + --- + ^» (i = 2,3,--- ,iV), 

2 ^ fL\ 

tpi = PiSi + Pi-iSi-i - j^^^ PiSi + ( ]y ~ ^M ^(^)- 

1=1 ^ ^ 

To determine ai, we add one more linear equation of the form 'Yl,i=i ouPi = P-i which is 
independent of the above A^ — 1 equations. Since 

N N 



i=l * i=2 * 



we obtain ai = {P — Q)/R. Next, we propose three candidates for each {pi} and P, and 
choose one of them in the following way: 
Candidate 1. We put 



N 



Pi 



Ti+i-Vi ,. ^ ^ ,n „ sr^ ri+i-2ri + ri_i 



^.^±V^ (^ = 1,2,---,A^), P = Y.Pl 



i=l 

and from (1551) we calculate R and Q. We denote this Rhy Ri. If the above equation 
holds, then err^i = and A = — X]j=i /^«*'"« hold in f l3ip . However, if distribution of grid 
points are almost uniform, then the above equation is almost nothing. Therefore we need 
another candidate. 

Candidate 2. For the z-th quantities Fj defined on Si and Gj defined on S*, we define 
the average along V such as 

N N 

{\=) = JY.^^n' (G)* = -5^G,r;. 

i=l i=l 

Since L = ^j^i^j = J2i=i''^h ^^ have (1) = (1)* = 1. Moreover, for a* = {ai + ai_i)/2 
defined on Si, the relation (a)* = (a*) holds. 

The second candidate of linear equation is the zero-average (a)* = 0, that is, pi = r* 
for i = 1,2,- ■ ■ ,N and P = 0. From this and fl33l) we calculate R and Q. We denote this 
R by i?2- 

The purpose of this section is to present numerical simulations of the geometric flow 
evolving according to the evolution equation ( !T5|l . Before we introduce the third candi- 
date, we calculate the discrete version of (ITSll . Let the total interfacial energy be 

N 

L^{t) = ^(j{vi{t))ri{t). 

i=l 

The time derivative of Tj = (cos z/j, sin Vi)^ is Tj = z>j ATj. Then we have rjZ>j = (ij — ij-i) ■ 
A^,-, from which it follows that 



N N N 

~iai, 



La = ^{(j'iUiri + aif-i) = - ^ kai/3*ri + ^ k^ip, 

i=l i=l i=l 

Pi = i.(^i + (^i+l)^i + i^i ~ (^i+l)Ci- 
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Here ai = a{ui), a[ = a'(i^j), and k^i is discrete version of the weighted curvature in the 
following sense: 



^i+i-'^Li , o-j+itj + cri(ti + tj_i) +crj_it 



i-l J. ^ ri 






"-"(TZ ^I'^l-j ^% 



2{ti + ti-i) 2{ti + ti_i) ' * 2' 



and 5j is discrete weight of 5{v) = o""(z/) + cr{i>) at u = z/j. Note that 6i — )■ 5(z/i) holds 
as (pi, (pi-i — )■ formally, and even the case where ti + ti_i = 0, k^i is well-defined, since 
ki = (ti + ti-i)/ri and then denominator of k^i is 2rj. 
We obtain 



A:S = ^ fi,^ _ ^A] , L. - %i = - y f A;.,, - ^) p*r, + err 
dt A A \ 2A J ' 2A ~t\ 2Ay ^* 

PiQi + ^ I 2^ A ^ 2^ "i^i — ^ 

Candidate 3. We put 

Ti+i-ri 2A. ,r-^ ri+i - 2ri + ri_i 
Pi = Si —Pi (« = 1,2, ■■■ ,iV), P = ^/3^ , 

^ ^^ i=i ^ 

and from (!33|) we calculate i? and Q- We denote this R by -R3. Let the i-th normal 
velocity defined on Si be 

I3* = k^,-^ (^ = l,2,---,iV), 

which is discrete version of flTSj) . If we choose the candidate 3 and its hold exactly, then 
err^atjo = holds and we obtain 

dLl 2L.^f, LA' 

Ti < 0. 






dt A A ^^\ "' 2A 

1=1 ^ ' 

Choice one from three candidates. We choose candidate number / satisfying 
\Ri\ = max{|i?i|, \R2\, l-RsD- 

Remark 3 //, for definition of the normal velocity defined on S* , we use 
P:n + /3*+ir.+i 



A 



2cir* 



instead of (E^j, then k^ri can not be divided into weighted part Wi and the curvature ki 
Moreover, if we use 



instead of (E^j; then k^i can be divided into weighted part and the curvature. However, 
we have err^ 7^ 0, that is A = — ^^^ 
ri = L/N holds for all i = 1,2,- ■■ ,N. 



we have err^ 7^ 0, that is A = —Yli=il^i^i ^^^■^ '^^^ ^^^^ even if uniform distribution 



On a gradient Bow of plane curves minimizing the anisoperimetric ratio 18 

Step 2: Discretization in time. We use semi-implicit scheme for discretization of 
( !29l) . Next we develop expression ( l29l) as follows: 

X, = aiT* + (3iN* 

1 foii f3i\ 1 fa, f3i\ 

Here we have used the relation 

' 2q ' ^^^ - ^^ 2 + ""^ 2 ~ 2si ■ 

Let /i be a parameter satisfying /i = if mini<i<Ar \si\ = 0, otherwise fi G (0, 1]. For the 
parameter /i G [0, 1] we use the linear interpolation of (IMl) and ( 1351) . since we can not use 
(ESI) if Si = 0. 

Put 6j = (1 — fi)si + fi/si for /i G (0, 1] and bi = s, for fi = 0. Then the evolution 
equation instead of f l29|) will be 

A. = ^^ - A&.) ^« + ^f + A&.) T,+i + Aq(1 - ^) ^-+^^-+\ (36) 

From Tj = (aij — a;j_i)/rj, we discretize (l36l) in time and obtain the following tridiagonal 
linear system under the periodic boundary condition: 

^ ^ = -a-^l^i + aoxr - a-,xi:i + ^M{1 - ^)^^^^^. 

Here oq = a_ + a+ and 



2rl \ d I 2rf 



J \ J 



i+l 



C: 



for z = 1, 2, ■ ■ ■ , A^ and j = 0, 1, 2, ■ ■ ■ . For the choice of /i we use 
min,- I Si I 

^^ = — 4^e[o,i]. 

maxj I s^ I 

Here and hereafter, miuj and maxj mean mini<j<7v and maxi<j<7v, respectively. Note that 
max,- Is;' I > holds for closed curves. 

In order to ensure solvability of the above linear system, we require a simple condition 
on the diagonal dominance. Adopting such a condition the adaptive time step Tj satisfies 

min,- r^ ,, s. 

r,- = ^^ ^ A > . 

' 2(1 + A)(max,|a^/c^|+ max, 1/3,^6^1) 

Simulation. In following all figures, for the prescribed r > 0, we plot every /^r discrete 
time step using discrete points representing the evolving curve. In every 3yur time step, 
we plot a polygonal curve connecting those points, where /i = [[T/r] /lOO] ([x] is the 
integer part of x), and T = 1.5 is the final computational time. We use u = 1000 as the 
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relaxation term in fl32|) . Note that if we use small u, then asymptotic speed for uniform 
distribution becomes slow, and for some choice of a area-decreasing phenomena do not 
hold near the initial time (cf. Figure H] (a)). 

Wulff shapes and area-decreasing phenomenon. If cr(i/) = 1 + e cos(rrai/) is the 
anisotropy density function of the degree m then, by using ([7]), we obtain the explicit 
expression for the mixed anisotropy function: 



a 



a/tt a + \/\W„\ = y/n 1 + 5 cos mu + \ 1 — e 



,m? — 1 



In order to verify the area-decreasing and the length-increasing phenomenon at the initial 
time as in Theorem HI we use dW^ as the initial curve. Figure [3] (a) indicates dW^ 
with m = Q. Its discretization is given by the uniform A^-division of the u-range [0, 1]. 
Figure E] (b) indicates the same Wulff shape, but the grid points are distributed uniformly. 
Figure [3] (c) indicates the time evolution starting from (b). 




Figure 3: (a) The Wulff shape dW^ with A^ 
and its time evolution (c). 



120 points, (b) its uniform parameterization 



Although deformation oi dW^ is very small (see Figure [3] (c)), the area-decreasing and 
the length-increasing phenomenon can numerically verified by using the aforementioned 
numerical discretization scheme. The behavior of the enclosed area and total length 
of curves evolved from the initial Wulff shape dW^ with the mixed anisotropy density 
function is shown in Figure HI 

Initial test curves. As initial test examples we use the boundary dW^ of the Wulff 
shape as well as the following initial curves x{u,0) = {xi{u),X2{u))^ {u G [0, 1]) parame- 
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Figure 4: Initial decrease of the enclosed area A and increase of the total length L. 



terized by 



(aj 
(b) 
(c) 
(d) 

(e) 



ellipse : X\{u) = acos27rn, X2{u) = 6sin27ru, 

dumbbell shape : Xi{u) = cosz, X2{u) = 2.0 sin z — 1.99 sin'^ 2, z = 27iu, 



Xi[Uj 



cosz, X2(m) = 0.7sin2 + sinxi 



X 



3' 



X3 = sin(32;) sins;, z = 2'itu, 



Xi{u) = 1.5 cos 2, X2{u) = 1.5(0.6 sin 2 + 0.5X3 + 0.4 sin 0:4 + 0.1 sin X5) 
Xs = sin(32;) sinz, X4 = 2x1, ^5 — 3e~^^, z = 27ru, 



thin-dumbbell shape : 

'xi(m) (0<m<0.5) 



XiiU 



Xi[U) 



—Xi{U - 
Xi{u) 

-xi(0.5 



0.5) (0.5 < M < 1 

(0 < M < 0.25) 
(0.25 <u< 0.5) 



, X2{u) 



X2{Uj 
-X2{u 



u 



X2[U) 



(0 <M<0.5) 
0.5) (0.5 <u<l) ' 

X2(m) (0 < m < 0.25) 

£2(0.5 -m) (0.25<m<0.5) 



- { \ - fbeam + rad{l + cos(8(7r - e^)u)) (0 < n < 0.125) 



X2[U) 



earn + rad{l - cosee))(l - Au) (0.125 < n < 0.25) 

radsm{8{7T-ee)u)) (0 < m < 0.125) 

e (0.125 <u< 0.25) ' 



where beam, > and < e < rad are parameters defined below and O,, = aicsin{e /rad). 

In all examples, the initial discretization is given by the uniform A^-division of the 
M-range [0, 1]. Figure [5] indicates numerical simulation with the initial curves. We choose 
several peaks of a such as (a) m = 2, (b) m = 3, (c) m = 4, (d) m = 5, and (e) m = 6. 

Breaking of comparison principle. In Figure [6] we plot the initial curve consisting 
of the union of the boundary of a Wulff shape dW^j (with a having degree m, = 3) touched 
by a circle with a sufficiently large circle with a radius r > 1 (see section 7). As soon as it 
evolved by the anisoperimetric ratio gradient flow it intersects the stationary Wulff shape 
dWa- Numerically computed examples displayed in Figure [6] (a) and (c) correspond to 
those of the conceptual Figure El 
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Figure 5: Evolution of curves starting from the initial curves with various choice of peak 

of <T. 
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Figure 6: Anisoperimetric ratio gradient flow breaking comparison principle. 



